FIRST-PASSAGE PERCOLATION ON WIDTH-TWO STRETCHES WITH 
EXPONENTIAL LINK WEIGHTS 
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Abstract. We consider the first-passage percolation problem on effectively one-dimen- 
sional graphs with vertex set {l...,n} X {0, 1} and translation-invariant edge-structure. For 
three of six non-trivial cases we obtain exact expressions for the asymptotic percolation 
rate x by solving certain recursive distributional equations and invoking results from er- 
godic theory to identify x as the expected asymptotic one-step growth of the first-passage 
time from (0, 0) to (n, 0). 



1. Introduction 

Let G = (V, E) be a graph with vertex set V = V(G) and unoriented edges E = E(G) c 
V 2 . Two vertices u, v e V are said to be adjacent, for which we write u ~ v if (m, v) e E 
and the edge e = (u, v) is said to join the vertices u and v. Assume there is a weight- 
function w : E — > R. For any two vertices u, v e V a path p joining u and v in G is 
a collection of vertices \u = po,P\, ■ ■ ■ ,p n -\,p n - v] sucn that p v and p v+1 are adjacent 
for < v < n; a path p is called simple if each vertex occurs in p at most once. To a 
given path p we associate the set of its comprising edges p = {(p v ,p v +i) : < v < n). 
The weight w(p) of a path p is then defined as Yieep w(e). We define dc : V x V — > R 
by dc(u, v) = inf {w(p) : p a path joining m and v in G), and call dc(u, v) the first-passage 
time between u and v. A path p joining u and v in G, such that w(p) = Jc( M > v ) is called a 
shortest path. Throughout this work we assume that G is finite and connected and we will 
be interested in the case that the weights w(e), e e E, are random variables; the goal is then 
to make probability statements about first-passage times or related quantities. In general 
we note here that a shortest path p = \u = po,...,p„ = v} is always simple and that each 
sub-path {pit, . . . ,pi), < k < I < n, is also a shortest path. Moreover, for continuously 
distributed, independent edge weights the shortest path between any two vertices is almost 
surely unique (Kesten, 1986). 

The typical first-passage percolation problem is based on two-dimensional regular in- 
finite graphs G with vertex sets V(G) = Z 2 . The edge weights w(e), e e E(G), are i.i.d 
random variables in L 1 with some common distribution P such that P-almost surely w(e) 
is positive. Let / m _>„, < m < n, denote the first-passage time from (m,Q) to (n, 0) subject 
to the condition that the contributing paths consist only of vertices with first coordinate v, 
m < v < n and write /„ := Iq^„. By Lemma 4.2, Zrj->* - fo->m + hn^n and the theory of 
sub-additive processes (Proposition 4. 1 ) implies that lim„^oo -l„ exists and is almost surely 
constant. This number, which depends only on the graph G and the distribution P, is de- 
noted by x(G, P) and called the (asymptotic) percolation rate or time constant. The explicit 
calculation of the percolation rate even for the simplest regular infinite graphs and distri- 
butions P is characterized in Graham et al. (1995, p. 1937) as a "hopelessly intractable" 
problem. 
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2. The model 

In this work we first focus on the first-passage percolation problem as described in 
the last paragraph for certain families of regular, effectively one-dimensional graphs with 
independent, random hypweights. Let V» = [0, 1, . . . , n} X {0, 1} and 

<V n ={((/, 0),(/+ 1,1)) :0<i<n] 

OV„ ={((/, 1),(/+ 1,0)) : < / < n) 

X„ ={((/, 0),(/+ 1,0)) : < i < n) 

y„={(0M),(j+l,l)):0< i<n] 

Zn ={((i, 0),(/,l)):0</<«). 

To each subset <5 c {'V, *W, X, J/, Z] corresponds a family of graphs 

Q & — {(?f}„(=N = {(Vji>£n)}neN 

with vertex sets V„ and edges S n = Uxee For eacn g ra P n the edge weights are indepen- 
dent exponentially distributed random variables which are labelled V,, Wi,X t , Y h Zj in the 
obvious way. By time-scaling it is no restriction of generality if we assume the parameter 
of the edge weight distributions to be unity. We denote the probability measure by P, that is 
dP(x) = e~ x dx. We point out that by construction, for any selection of edges S, G% is a sub- 
graph of Gf, ,, an observation which forms the basis for our inductive argument described 
below. The method is based upon Flaxman et al. (2011) where it has successfully been 
employed to compute the asymptotic percolation rate on where the edge weights 

are taken to be independently one with probability p and zero with probability 1 - p. They 
also consider continuous edge-weight distributions and show that one can replace those 
by suitably chosen discrete ones to obtain arbitrarily good, yet approximative values for 
the time constant. We do not adopt this method here but rather work explicitly with the 
continuous distributions. 

3. Results 

It turns out that for only six families of graphs Q E the first-passage percolation is non- 
trivial. For three of these six, namely for Q^^'^, gWW.XJf] an( j gl'WjCj/.Z) we (j er i ve 
expressions for^ which seem not to have been known so far (see Table 1). We are con- 
fident that our method works for the three remaining families, gWWW 7 gWfWji^ an( j 
g{ r v, r w,x,y,Z] ^ as we u Increased complexity in the calculations involved, however, prevents 
us from obtaining analytic results; the numerical values of x f° r theses cases, which are 
also recorded in Table 1, were obtained by means of simulations. 

4. S UBADDITI VITY AND ERGODIC THEOREMS 

In this section we present some theoretical results about first-passage times. In particular 
we show that the asymptotic percolation rate lim^oo -l n on graphs with vertex set V„ := 
{0, ...,«} X {0, 1} and translation-invariant edge structure is, with probability one, equal to 
a deterministic number and provide a formula for its computation (Theorem 4.4). 

Proposition 4.1. Let G be the graph G% for some & c { < V, r W,X,}f,Z}- If the edge 
weights w(e), e € &„, are independent, identically ^-distributed, positive and integrable 
random variables then almost surely lim„_, M -/„ exists and is equal to a deterministic num- 
ber x(G,V). 

To prove this we need a lemma about a property of first-passage times, called subaddi- 
tivity. 

Lemma 4.2. For any weighted graph G, the function dc ■ V(G) 2 — » R is subadditive, that 
is dc(u, w) < dc(u, v) + dc(v, w) for all u,v,w G V(G). If in particular V(G) = V n , this says 
lo-> n - lo->m + lm->n holds for any < m < n. 
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Proof. Subadditivity of d c is a consequence of the simple observation that the set of all 
paths joining u and w in G is a superset of the set of all paths joining u and w in G and 
containing v and that one can combine any two paths p M_>v joining u and v and //'^" joining 
v and w to get a path p"^ w joining u and w that satisfies w (p u ~* w ) = w (p M_>v ) + w (p v_>w ). 
The remark about the special case G = Gf is clear. (Take m = (0, 0), v = (m, 0), w = 
(n,0)) □ 

To conclude from this subadditivity property that the asymptotic percolation rate is al- 
most surely a deterministic number we need the following subadditive ergodic theorem 
due to Liggett, who generalized a result of Kingman. (See Durrett (1991, theorem 6.1) for 
a proof.) 

Theorem 4.3 (Liggett). Suppose a family of random variables X - {X„ h „ : < m < n) 
satisfies 

(i) X^ n < X() 

(ii) For each k € N, {X n k^ n +i)k)n>x> is a stationary sequence. 
(Hi) The distribution of (X m m +k)k>o does not depend on m. 
(iv) E [Xq J < oo and for each n e N, E [Xo ; „] > yo holds with yo > -oo. 

Then X — lim,,-^ Xo jfl /n exists almost surely. 

In the formulation of (iv) we used, as we will often do in the following, the notation x + 
as a short-hand for max{0, x}. We can now give a proof of the assertion that the asymptotic 
percolation rate almost surely equals a non-random number. 

Proof of Proposition 4.1. We will first argue that the family \X mn = Z m _,„ : < m < n] sat- 
isfies the conditions of Theorem 4.3. In Lemma 4.2 we have seen that (i) holds, (ii) and (iii) 
follow from the translational invariance of the graph and the fact that {w(e),e e E} is i.i.d. 
The last condition, (iv), holds by assumption and in particular for the exponential distribu- 
tion. We can thus conclude that l n /n — > x almost surely as n — » oo. To see that^ is indeed 

constant, we enumerate the edges \e e E} in some way, say e\ , ei, Since^ is defined as 

lim, MOO l„/n and each /„ is a sum of only a finite number of edge-weights, namely of a sub- 
set of those contained in Gf, x is measurable with respect to Hn>i o"(w(e K ), w(e„+i), . . .), 
the tail-cr-algebra of the i.i.d. sequence w(e\), w(e2), Kolmogorov's Zero-One law (see 
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for example Breiman (1968, Theorem 3.12)) then implies that for any Borel set B the prob- 
ability that x takes a value in B is either or 1 which means that x is almost surely equal 
to a deterministic number. □ 

In view of our next theorem, which provides an explicit formula for the percolation rate 
X, we could have done without Proposition 4.1 and without invoking Kingman's general 
subadditivity result. We chose to include it, however, in order to distinguish this universal 
aspect of first-passage percolation theory from properties specific to our model. The result 
is the following: 

Theorem 4.4. Let A„ :— l„ — l„-\ and assume there exists an S -valued ergodic Markov 
chain (M n )„>\ and a measurable function f : S — * R satisfying A„ = f(M n )for all positive 
integers n. Then, almost surely, 

(1) X = lim - = f f(sMs), 

n-wo n J s 

where n is the unique invariant distribution of (M„) n >\. Put differently, x — E[A], where 

D 

A is the weak limit of the sequence (A„)„>i, i.e. A„ — > A. 

Remark 4.5. It will be established in Lemma 5.6 that there indeed exists a Markov chain 
(M„)„>i and a function / satisfying the conditions of the theorem. 

Proof. As an instantaneous function of the ergodic Markov chain (M„)„>\ the sequence 
(A„)„>i is ergodic as well. It follows that the time average i Y!l=\ Ay converges to E [A]. 
Since clearly l n = £" =1 A„ the claim follows. □ 



5. Calculations for Q^mOi 

In this section we prove our results about the time constants for first-passage percolation 
on width-two stretches (Table 1) in the case of S — \X,}f,X}- The calculations for the 
other two cases mentioned are completely analogous and can be found in Schlemm (2008). 
We denote the length of the shortest path from (0, 0) to (n, 0) by l„, the length of the shortest 
path from (0, 0) to (n, 1) by l' n and we let A„ = l' n — l n . Our first goal is to find a recurrence 
relation between the distributions of A„ and A„_i. 

Proposition 5.1. The sequence of random variables (A„)„>o is a real-valued Markov chain 
with initial distribution P (Ao < d) = \ — e~ d and transition kernel 

1 d < A 5 < d, 

e -(s-d) d < A 5 > d, 

e-M d = 0, 

e - {d - S) d>QA6<d, 

1 d > A 5 > d. 



(2) K(S,d) = e 



-\d\ 



Proof. It is clear that we must try to express p„, the probability density function of A„, in 
terms of p„_i, the density function of A„_i. To achieve this we define events on which A„ 
is a deterministic function of X n , Y n , Z„ and A„_i, which amounts to separately taking into 
account the different possible behaviours of the last step of the shortest path. Consider first 
the shortest path p ^ 1 from (0,0) to (n, 1), the last edge of which must be either Y n or Z„. 
Similarly, the last edge of p°^ l> , the shortest path from (0, 0) to (n, 0), must be either X n or 
Z„. We note that almost surely, the edge Z„ is not part of both p ^ 1 and p 0_>0 because this 
would contradict the uniqueness property mentioned in the introduction, so that we have 
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to consider three cases: 
(3a) 



(3b) 



- [Y n + 4-1 - Z n + X„ + l n -\ A X„ + /„_i < Z n + Y„ + l' n _ l 
= {Z n >\X n -Y n -A n . l \] 

= {z„ + x n + /„_] < f„ + 4_i a x„ + i n -\ <z„ + y„ + i' n _ x J 

= {F„>Z„+Z n -A n _ 1 } 



(3c) — = {y„ + /;_] < z„ + x„ + a z„ + y„ + < x n + j^) 
= {X n >Y n +z n + A n . l ] 

In particular we have the relation 

(4) A n = min{A„_i + Y„, X n + Z n ) - mm{X n , An—I + Y n + Z n ). 

The reason why the pictographs above have been chosen to represent the different events 
is the following: solid lines correspond to edges being used by the shortest path to either 
(«,0) or (n, 1), double solid lines represent edges being part of both these paths, while 
dashed lines stand for unused edges. Using these events we now compute the cumulative 
distribution function of A„ as 

(5) P(A„ < d) = P({A„ <4nj) + P({A„ <d}n"3) + P({A„ < d} nZJ) . 

On the event l, A„ is equal to A„_i + Y„ - X„, so the first term is 



P(A„ <dn^) = J ddpn-M c/F 3 (^:,y, Z) li^l.v-v-tfKlitf+v-.vsrf}- 
Rewriting the indicator functions as integration bounds we obtain 

dSp n -m d¥(y) dP(x) dP(z) 

JO J(t5+v-rf)+ J\x-y-S\ 

and after doing the tedious but easy x-, y- and z-integrals we arrive at 



>(a„ <dn. )= f d6p n -i(S)Gi(6,d) 
Jr 



with 



Gi(M) 



1 



e s (l + 2(d-6)) d<QA6<d, 

e 2(d - 6) d<0A6>d, 

e s {2 -26- e- 2d ) d > A 6 < 0, 

e- & (2 + 25 - e - 2{d - 5) ) d>Q AQ <6 <d, 

e~ s (l +26) d>0A6>d. 

One can also confirm this computation and similar ones which follow by use of a computer 
algebra system. In fact, we used Mathematica to verify our results. On the event 3, A„ is 
given by Z„, so for the second term in (5) we obtain: 

P(A„ <dfO= f d6p n -i(6) f d¥\x,y,z)I tv >, +z - di l {z < d } 
Jr Jr 3 

• y-»00 y-»00 

d6p„^{6) I dP(x) ( dP(z) I dP(y) 

l Jo Jo J(x+z-6)+ 

d6p n -i(6)G 2 (6,d), 
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where the function G 2 is given by 

e s {\-e- 2d ) 6<0, 



1 

G 2 (6,d) = -I|(/>o) 



4 - e~ s (3 + e - 2(d - s) + 26) 0<6<d, 
4 - 4e- d - 2de- s 6 > d. 



Finally, on the event A„ is equal to -Z„, so we compute for the third term in (5) 

P(A„ < d n H) =Jddp n - 1 (d) ^ dP 3 (x,y,z)l u >y + z + s)h<-d) 

d6p-M5) ( <ff(y) ( dP(z) I c/P(x) 

Jo J(-<0+ J(.v+z+<5) + 



d8p n -\($)G-$(6,d) 



and 



G 3 (<?,</) 



4e* + e' 5 (-2(t/ - 5) - 3) d<0A6<d, 
1 j e- 2 ^-* d < A 6 > d, 

4 |4 + e <5 (2£-3) d > A 6 < 0, 

e~ s d > A 6 > 0. 

It is interesting to note that, due to the symmetry of the events (3b) and (3c), the sum of 
G2(6, d) and G 3 {-6, -d) does not depend on d; explicitly it holds that 

1 Wl fl 6<0 
G 2 (6, d) + G 3 (-6, -d) = -e~ lsl \ . 

4 ]4e s -26-3 6>Q 

Putting everything together it follows that 

P(A„<d) = I d6p„- l (6)[G l (6,d) + G 2 (6,d) + G 3 (6,d)]. 
Jr 

Differentiating this equation with respect to d and interchanging differentiation and inte- 
gration on the right hand side we obtain 

p n {8)= f d6p„- 1 (6)dAG l [6,d]+G2(6,d) + G 3 (6,d)], 



so the transition kernel K : R -> R is given by K(6, d) = d d [G x [6, d] + G 2 (6, d) + G 3 (6, d)] 
and the claim follows from basic computations. □ 

From the explicit description of the transition kernel K{6, d) we can deduce useful prop- 
erties of (A„)„> : 

Lemma 5.2. The Markov chain (A„)„>o is ergodic. In particular, as n — > oo, (A„)„>o 
converges in distribution to a non-degenerate limiting random variable A with probability 
density function p^ — lim^oo p n . 

Proof. It is enough to note that the kernel function K from Proposition 5.1 is strictly 
positive and continuous. The claim therefore follows from an extension of the Perron- 
Frobenius theorem to continuous transfer operators (see lentzsch, 1912). □ 

Corollary 5.3. Let p M = lim^oo p„. Then 
(i) The density p^ of the stationary distribution o/(A„)„>o satisfies the integral equation 



(6a) 



Xi r*oo 
d6 Poo (6) + e ld d6 Poo (6)e- 5 d<0 
OO xJd 

Xd poo 
d6 Poo (6)e 5 + e- d d6 Poo (6) d>0. 
00 %Jd 
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(U) The density p^ is an even function, that is Poo(d) = p m (—d) for all d e R. 

Proof. For the first claim we observe that p^, satisfies the integral equation p^id) = 
J K d6 poo(6)K(6, d) with the kernel K given in (2). The second assertion follows from the 
symmetry K(5, d) = K(-5, -d). a 

In order to solve this integral equation we transform it into a differential equation. 

Lemma 5.4. The density p M of the stationary distribution o/(A„)„>o satisfies the differen- 
tial equation 

(7a) p^(d)=-[2 + e d ]p oa (d)+3p' co (d) d < 0, 

(7b) p'M=- [2 + e- d ] pM) ~ 3p'Jd) d>0. 

Proof. We only prove the claim for d < 0. The case d > can be shown in the same way 
or one uses (ii) of Corollary 5.3. In the case d < we know that p M is a solution to the 
integral equation (6a). Differentiating this equation with respect to d and using (6a) again 
we obtain 



p'^id) =e d pUd) + e d J d6p 

%J — CO 



(8) 



Jr>oa 
d5 Poo (5)e- s 
d 



-2p oa (d)-e d | ddpoo(S). 

Differentiating again and eliminating the remaining integral via (8) yields 

r d 

pl(d) ^2p' m (d) - e d pUd) - e d dSp x {S) 

%J —CO 

=3p'Jd)-[2 + e d ]p so (d). 

□ 

Proposition 5.5. 77ie density p M of the stationary distribution of the Markov chain (A„)„>o 
is 

(9) p„(d) = -^—e-^h(2e-^), 
where J v are Bessel functions of the first kind. 

Proof. For d < 0, we write p x as p 00 (ii) = e^p™ (2ei £/ ) with a function p x which is to 
be determined. If this ansatz is inserted into (7a) and if we replace d by 2 log | we obtain 

z 2 ^£r(z) + z^iz) + (z 2 - l)poo(z) = 0. This is the differential equation the solution to 
which are by definition the Bessel functions of the first and second kind, denoted by J v 
and Y v , respectively (see Abramowitz and Stegun, 1965, equation 9.1.1). In this particular 
case v = 1. One can apply an analogous method for the case d > to find that the general 
solution of equation (7) is given by 

(10) Poo(d) = e-^ [c, J, (2e-'^) + c 2 Y 1 (2e~^)] . 

In order to determine the two constants c\ and c^ we insert this general solution into the 
integral equation (6) and take the limit d — > (it does not make a difference whether we 
take d < or d > 0). We choose d < and obtain 



X0 poo 
£%«(£)+ J ddpccid^ 6 , 
oo Jo 
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in which we replace pco(£) by the expression given in (10) to get 



c 1 J I (2) + c 2 Y 1 (2) = Cl 

+ c 2 

=d [J 2 (2) + J (2)]+ C2 



X0 poo 
dSe^- d h{2-e } i d ')+ I dSe^- d h {le'^e 
CO J() 

X0 r*co 
ddei d Y 1 (2ei d )+ d6e^ d Yj {2e^ d ) 
co %JQ 



1 



+ Y 2 (2) + -- +Y (2) 



Since Ji (2) = J (2) + J 2 (2), Yi (2) = Y (2) + Y 2 (2) and -l/n + it follows that c 2 = 0. 
The requirement that p x integrate to one implies 



= f d6e-^h(2e-n=\ f dyy 2 h(y) = 2J 2 
Jr. z Jo 



(2). 



Lemma 5.6. There exists anRx R^-va/wec/ ergodic Markov chain (M„)„>i one/ a function 



f 



. such that A„ = f(M n )for all n > 1. 



Proof. Set M„ = (A„_i,X„, T„,Z„) and f(6, x,y,z) — min{x, 6 + y + z). It is then clear that 
f(M„) = A„ (see (12)) and (M„)„>i is Markov because A„_i can be written as min{A„_ 2 + 
F„_i , X„_! + Z„_i } - min{X n _! , A„_ 2 + F„_! + Z„_i }. (cf. equation (4) in the proof of Propo- 
sition 5.1.) Ergodicity is a direct consequence of Lemma 5.2. □ 

We will now use this result to compute the distribution of A := lim,,-™ A„, the expected 
value of which is, by Theorem 4.4, the percolation rate we are looking for. 

Lemma 5.7. For each natural number n, the density r\ n of the distribution of A„ is tj„(1) — 
J R d6 p„-x(8)Q(6, 1) with Q : R 2 -> R given by 



(11) 



Q{5,l) = e- 



e\l -5) I < A 5 < I, 

e -«-<5)(l +2(1-6)) 1>0A6<1, 

1 I > A 6 > I, 

otherwise. 



In particular, ^(1) = J R dS p oa {6)Q(8,l). 

Proof. As before we define events for each possible behaviour of the last step of the short- 
est path from (0, 0) to (0, n): 



■ \X n + l n -\ < Z n + Y n + t_j = \X„ < A„_i + Y„ + 

■ [x n + /„_, > z„ + y„ + = {X n > A n _, + Y„ + Z n ) 



Clearly, 
(12) 



A„ = min{X„, A„_i + Y„ + Z n }, 



i.e. on the first of the two events, A„ is given by X„ while on the second it equals Y„ + 
Z„ + A„_i. Now the procedure continues very similar to the proof of Proposition 5.1. We 
compute the cumulative distribution function of A as 



(13) 



'(A„ < o = p({A„ < /( n ) + p({A„ </)nJ), 
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because on each of these two events we have an explicit expression for A„. The first term 
leads to 



t ({A n <l}nZ)= f dS Pn -m f dF\x,y,z) 
v Jr Jr 3 

dSp^m d¥(z) dP(y) 

Jo Jo Jo 

d5p„- X (6)Px(6J), 



{x<d+y+z)Hx<l} 

min|5+v+z,/) + 



to 



where the function P\ is the result of the x-, y- and z-integration and given explicitly by 

1 e -2i+s p/ (3 _ 25) _ 3 _ 2(/ _ 6 )j 5 < o, 



Pi (5,0 =I{/>o t 



[4 - e" 5 - e- 2l+s (2(l - 6) + 3)] < 5 < /, 
1 - e"' 6 > I. 



Very similarly we obtain for the second term in (13) 
P({A„ <l}n ll) = J^ d6p n -i(S) J~^dF 3 (x 

. )>, Z) I{ x >d+y+z)hs+y+z<D 



Jo Jo 
ddp n - 1 (S)P 2 (S,[), 



dP(y) 

J(<5+y+z) + 



to 



with P 2 given by 



1 -e'-'\l -6 + 1) 1<0a6<1, 

1 - \e & [3 - 26 + e- 2/ (l + 2(/ - 5))] Z > A 6 < 0, 

i g -2/-a yii _ e 26 (l + 2Q _ £))] />0a0<<5</, 

otherwise. 



Putting things together, it follows that P ( A„ <l) = f R d6p ni (6) [Pi (6, I) + P 2 (6,l)}. Taking 
the derivative with respect to I and interchanging differentiation and integration on the 
right hand side we obtain Q{6, I) = d[ [P\(6, I) + P 2 (6, 01- The statement about the relation 
between the stationary densities 77^ and p x is clear. □ 



Theorem 5.8. The percolation rate x on $ x #& £y [| - Jj^U a* 0.68 .... 

Proof. We have already argued that the time constant is the expectation of A, that is 
X = L dllr/coil). From Lemma 5.7 we know the explicit form of r/ca so we obtain^ 
I . d6p oa (6) L dllQ(6,l), where the change of the order of integration is easily justified 
using Fubini's theorem. A direct calculation shows 



dllQ(6,l) = 



1 8 + 46 - e s (5 - 26) 6<Q 



4 4- 



6> 



and the percolation rate is 
1 



X = 



2J 2 (2) 



\ dSe* 5 h{^ s ) [8 + 46-e s (5-2d)] 

Jr_ 



f d6e-^J l (2e-^)U-e- s ] 

Jr + 



1 

: 2J 2 (2) 



4J 1 (2)-7J (2) 4J 2 (2)-J (2) 



Ji(2) 
2J 2 (2)' 
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6. Discussion 

The differences between the omitted computations for the cases {*V, *W, X, J/} and 
{ r W, X, if , X) and the ones presented are only a matter of degree, not of kind; for instance, 
there are not only three different cases to consider when proving the analogues of Propo- 
sition 5.1, but rather eight and six, respectively. Moreover, for the case { < W, X, if, X) 
the stationary density p x of (A„) n >i (cf. Proposition 5.5) is not symmetric. It is also 
the determination of this stationary distribution where things get more difficult with the 
unsolved cases: the characterizing differential equation (cf. Lemma 5.4) then becomes 
non-local, more specifically it involves both Poaid) and p^i-d), as well as their deriva- 
tives, at the same time and so one can not solve it separately for d > and d < as we 
did. It might be possible to remedy this by computing p^ not as an eigen function to the 
transition kernel K itself but rather as an eigen function to its second convolution power 
K (T> (5, d) — C dcrK(6, cr)K(cr, d). As the method used in this paper is very similar to that in 
its antecedent Flaxman et al. (201 1) its range of applicability is also essentially the same. It 
would be a natural generalization to consider graphs with vertex sets {1, . . . ,«} X {0, . . . , k] 
for some integer k > 1 and for the directed first-passage percolation problem, similar 
techniques as used here apply equally well to this more general set-up. However, the com- 
binatorial difficulties arising from the need to explicitly keep track of the shortest paths 
on G„+i\G„ seem very hard to overcome and for undirected percolation a similarly easy 
recursive argument as in the case k = 2 is not an option. 
Many thanks go to Balint Virag for very helpful advice. 
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